Can vacancies lubricate dislocation motion in aluminum? 
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The interaction of vacancy with dislocations in Al is studied using the Semidiscrete Variational 
Peierls-Nabarro model with ab initio determined 7-surface. For the first time, we confirm theo- 
retically the so-called vacancy lubrication effect on dislocation motion in Al, a discovery that can 
settle a long-standing controversy in dislocation theory for fee metals. We provide insights on the 
lubrication effect by exploring the connection between dislocation mobility and its core width. We 
predict an increased dislocation splitting in the presence of vacancy. We find that on average there is 
a weak repulsion between vacancies and dislocations which is independent of dislocation character. 



Defects and their mutual interactions dominate the 
properties of materials that host them. Vacancies as 
point defects, have long been known to strongly interact 
with dislocations (line defects), and the study of their in- 
teractions represents one of most challenging problems in 
material science and engineering Q . More than a decade 
ago, Benoit et al. discovered an interesting phenomenon 
in ultra-high-purity aluminum deformed at low tempera- 
ture of 4.2 K. They observed a marked decrease of elastic 
modulus in cold-worked Al, which they attributed to va- 
cancy enhanced dislocation mobility in Al, a novel phe- 
nomenon they termed dislocation lubrication effect j^, [| . 
Corroborative experimental evidence led these authors 
to conclude that vacancies, generated from cold- work or 
irradiation, are solely responsible for the enhanced dis- 
location mobility. This indeed is a quite intriguing re- 
sult because traditionally vacancies are thought to lock 
dislocation motion by forming atmospheres around the 
dislocation [Q. Furthermore the lubrication effect may 
hold the key to resolve the long-standing controversy 
for Peierls stress (o~ p ) estimated from plastic deformation 
and from internal friction measurements. It is generally 
believed that kink pair formation (KPF) of dislocations 
is responsible for the Bordoni peaks observed in internal 
friction measurements of fee metals ||. However o~ p 
derived from KPF mechanism is in the order of 10~ 3 /i 
(/i is the shear modulus) far greater than the critical re- 
solved shear stress (CRSS) estimated from plastic de- 
formation experiments, which is around to 10~ 4 /i 
. This has been regarded as a serious problem because 
it casts doubt on the well accepted KPF theory for the 
Bordoni peak in fee metals. The controversy is particu- 
larly troublesome in light of the good agreement found in 
bec and ionic crystals regarding o~ p measured from inter- 
nal friction and low temperature CRSS experiments 
Therefore the vacancy lubrication effect may settle the 
controversy by proposing that vacancies strongly inter- 
act with dislocations and as a consequence lower their a p 
to the level that is consistent with the low temperature 
CRSS experiments and thus bridge the gap in a p 0. 
Finally if the lubrication mechanism turns out to be gen- 
eral, it may lead to an innovation in molding technology 



for materials with high a p by introducing vacancies to 
the materials Q . However interesting the lubrication ef- 
fect may seem to be, it has not been widely accepted, 
which in our opinion, is due to poor understanding of 
the phenomenon. In fact, there is no complete theoreti- 
cal work ever published to address this problem, to the 
best of our knowledge. Therefore it is the purpose of 
this paper to present the first ab initio study of the prob- 
lem. As we will show in the following, our calculations 
not only provide theoretical support for the lubrication 
effect, they also reveal other important difference in dis- 
location properties that are associated with the presence 
of vacancies. 

In this paper, we employed the recently developed 
Semidiscrete Variational Peierls-Nabarro (SVPN) model 
[E^lIJ in conjunction with ab initio determined 7-surfaces 
[11|. The SVPN model provides an ideal framework for 
multiscale simulations of dislocation properties, and it 
combines an atomistic (ab initio) treatment of the inter- 
actions across the slip plane and an elastic treatment of 
the continua on either side of the slip plane. The model 
has been shown to be quite successful in predicting dis- 
location properties by comparing its predictions against 
the direct atomistic simulations results [0, [nj. For ex- 
ample, by using the 7-surface calculated from Embedded 
Atom Method (EAM), we obtained a p for dislocations in 
Al that are in excellent agreement with that from direct 
atomistic simulations employing the same EAM potential 
fiof . More remarkably, a good agreement is also achieved 
for dislocations in Si || where the classic Peierls-Nabarro 
model fails owing to the its insufficiency to deal with nar- 
row dislocations, such as dislocations in Si. Since one 
has no a priori knowledge regarding the size of dislo- 
cations in the presence of vacancies, the SVPN model 
seems to be particularly useful to explore the interaction 
of vacancy with dislocations. Besides o~ p , the model can 
also provide reliable results for other dislocation proper- 
ties, such as partial separation distance and core width 
p0| , as they are compared to direct atomistic simula- 
tions |13, 14 . Thus the strength of this approach, when 



combined with ab initio calculations for 7-surface is that 
it produces essentially an atomistic simulation for dis- 
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location properties without suffering from uncertainties 
associated with empirical potentials. 

In the SVPN approach, the equilibrium structure of 
a dislocation is obtained by minimizing the dislocation 
energy functional |, Q 

Udisl — U e l as tic + V m i s fit + V 'stress + Kb 2 \llL, (1) 

where 

(2) 

U mis fit = J~] Axj(fj), (3) 

^. = -£^^=MM>, (4) 

i,i 

with respect to the dislocation Burgers vector density pi. 
Here, p^ 1 ' , p- 2 ^ and p\ 3 ' are the edge, vertical and screw 
components of the general Burgers vector density defined 
at the ith nodal point as pi = (fi — fi-i)/(xi — 
where fi and X i SX6 the displacement vector and the co- 
ordinate of the zth nodal point (atomic row). j(fi) is the 
7-surface that is determined from ab initio calculations. 

is the external stress component interacting with the 
corresponding Burgers vector density p"\ \ij is the dis- 
cretized elastic energy kernel, and K, K e and K s are the 
pre-logarithmic elastic energy factors H [l(J. The quan- 
tity L entering the last term is the outer cutoff radius for 
the configuration-independent part of the elastic energy 
p2| . We identify the dislocation configuration- dependent 
part of the elastic energy and the misfit energy as core 
energy, i.e., U core = U e i a stic + U mis fit- The response of 
a dislocation to an external stress is achieved by mini- 
mization of the energy functional at the given value of 
the applied stress. An instability is reached when an op- 
timal solution for the Burgers vector density distribution 
no longer exists, which is manifested numerically by the 
failure of the minimization procedure to convergence. o v 
is then identified as the critical value of the applied stress 
giving rise to this instability. 

In order to examine how vacancies change dislocation 
core structure by modifying atomic bonding across the 
slip plane, we carry out ab inito calculations for the 7- 
surface of Al with vacancies at the slip plane. Specif- 
ically, we select a supercell containing six Al layers in 
[111] direction with four atoms per layer, and remove 
one Al atom from the top layer (right below the desig- 
nated slip plane) of the supercell to simulate a vacancy 
concentration at 4 at.%. We should emphasize that 4 
at.% represents the vacancy concentration at the dislo- 
cation core region that we are interested, therefore it is 
much greater than the average vacancy concentration of 



the bulk material. The ab initio calculations are based on 
the pseudopotential plane-wave method with local den- 
sity approximation |l5j to the exchange-correlation func- 
tional |l(|. A kinetic energy cutoff of 12 Ry for the plane- 
wave basis is used and a fc-point mesh consisting of (8,8,4) 
divisions along the reciprocal lattice vectors is sampled 
for the Brillouin zone integration. Atomic relaxation is 
performed before we initiate the sliding. During the slid- 
ing process, atoms are allowed to move only along [111] 
direction while the atoms at the innermost two layers 
are held fixed. Volume relaxation is also performed for 
each sliding distance to minimize the tensile stress on the 
supercell. 

The ab initio determined 7-surface for Al with and 
without vacancies is presented in Fig. 1. In order to 
highlight the vacancy effect on 7-surface, we also sum- 
marize in Table I some important stacking fault ener- 
gies for both Al and Al+V (Al with vacancies) systems. 
These special stacking faults correspond to the various 
extremes along [121] and [101] directions of the 7-surface. 
It is clear that the vacancy lowers the intrinsic and un- 
stable stacking fault energy along [121] direction while 
increases the run-on stacking fault energy and unstable 
stacking fault energy along [101] direction. Therefore it is 
not immediately clear how dislocation core structure can 
be changed by vacancies, and a detailed analysis based 
on SVPN model is needed. Since the experiments Q 
have concluded that the change in elastic constants due 
to vacancies is not responsible for the observed lubrica- 
tion effect, we will simply use the experimental elastic 
constants of pure Al in our calculations for Al+V system 

Having determined all the necessary parameters enter- 
ing the model, we can study the interaction of vacancy 
and dislocations using the SVPN model. We select four, 
namely, screw (0°), 30°, 60° and edge (90°) dislocations, 
all with the same Burgers vector, but different orienta- 
tions, in our calculations. o v calculated from the SVPN 
model for the dislocations with and without vacancies 
is listed in Table II. a p represents the intrinsic mobility 
of a straight dislocation, and it relates to the kink pair 
formation energy, 2Wk, by 

where b is the Burgers vector and a the distance between 
neighboring Peierls valleys [|[ ^, |l7| . The correspond- 
ing values for Al are p = 28.8 GPa, b = 2.85 A and 
a = 2.47A. By measuring 2Wk (the activation energy for 
the Bordoni peak) in internal friction experiments, one 
can derive o v according to Eq. (5) for the relevant dislo- 
cation. For example, using the experimentally measured 
value of 2Wk (0.21 eV) for a screw dislocation, one ob- 
tains dp = 8xlO~V (224 MPa) for the screw dislocation 
in pure Al || [TtJ . This value of a v is in excellent agree- 
ment with our model result, 8.82xl0 _3 /z for the same 
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dislocation (Table II). Furthermore, the less definite mea- 
surement for the subsidiary peak (Bl peak) yields an ac- 
tivation energy ranging from 0.12 to 0.16 eV, which corre- 
sponds to a p in the range of 2.8 to 4.6x10 _3 /j, (80 to 130 
MPa) for the 60° dislocation. This value also agrees well 
with our result for the same dislocation (3.40x 10~ 3 ^i). 
The overall consistency between the theoretical and ex- 
perimental values for a p indicates the reliability of our 
model and establishes the basis for further study of va- 
cancy effect. 

When vacancies are introduced at the slip plane but 
are not adsorbed by a dislocation line, we find that o~ p 
for various dislocation is lowered by more than one or- 
der of magnitude (except for the edge dislocation), as 
shown in Table II. Therefore we have confirmed the va- 
cancy lubrication effect theoretically for the first time 
since its proposal. The fact that this lubrication effect is 
observed for various dislocations suggests a generic na- 
ture of the underlying mechanism. In order to shed light 
on this general mechanism, we have calculated disloca- 
tion core width which is defined as the atomic spacing 
over which the relative displacement of the dislocation 
changes from 1/46 to 3/4& fl(J. It is generally believed 
that o~ p is exponentially lowered with the increase of dislo- 
cation half- width according to the Peierls-Nabarro model 
0, H). The calculated dislocation half core width (£) 
is presented in Table II. It is found that in the pres- 
ence of vacancy, dislocation becomes 60% to 90% wider, 
which we believe is due to the reduced slope of the 7- 
surface along [121] direction as vacancies are introduced. 
Since lattice restoring force, represented by the slope of 
7-surface, is weakened by the vacancies, the repulsive 
elastic force resulting from the continuous distribution of 
infinitesimal dislocations dominates, leading to a wider 
dislocation core and therefore enhanced dislocation mo- 
bility. Although the vacancy lubrication effect may be 
qualitatively understood from above argument by a care- 
ful inspection of the 7-surface, one has to resort to SVPN 
model to obtain reliable values of a p in order to make a 
quantitative comparison. As shown in Table II, vacan- 
cies can bring a p down to the values derived from the 
plastic deformation experiments (10~ 5 to 10 -4 /i), there- 
fore bridge the gap for a p between the internal friction 
measurement and the measurement of CRSS at low tem- 
perature. 

In order to gain more insights into the interaction of 
vacancy and dislocations, we have calculated dislocation 
Burgers vector density for pure Al and Al+V, shown in 
Fig. 2. It is found that dislocations tend to dissociate 
more into partials in the presence of vacancy. This behav- 
ior is obviously associated with the fact that the intrinsic 
stacking fault energy is reduced by vacancy. The result 
cautions us to be more careful in interpreting TEM data 
for partial separation distances because accidentally in- 
troduced vacancies could change the result significantly. 
We have also calculated binding energies of vacancies to 



dislocation cores, summarized in Table II. The binding 
energy is defined as difference between dislocation core 
energy with and without the presence of vacancy. Over- 
all we find that the binding energies are not sensitive 
to dislocation character, and more importantly they are 
all marginally positive. The positive binding energy in- 
dicates that dislocation is energetically less stable with 
vacancies nearby. Our result for positive value of binding 
energy qualitatively agrees with the findings from atom- 
istic simulations carried out for Cu B. The small value 
of binding energy we found in Al may relate to the fact 
that the dislocations are dissociated into partials in Al. 
It is believed that the vacancy binding energy to dissoci- 
ated dislocations is considerably lower than that to the 
same dislocations in the non-dissociated condition |E1 . Of 
course one has to bear in mind that the binding energy 
we obtained represents an average value over all possible 
sites for a vacancy at the core of the dislocation. 

To conclude, we have studied the interaction of va- 
cancy with dislocations in Al using the SVPN model with 
ab initio determined 7-surface. We confirm the experi- 
mental finding of vacancy lubrication effect in Al. We 
propose that vacancies can weaken the lattice restoring 
force across the slip plane, which leads to a wider dis- 
location spreading, and thus higher dislocation mobility. 
We find that a p of the dislocations is lowered by more 
than one order of magnitude in the presence of vacancy, 
which bridges the gap between o~ p observed from different 
experiments, resolving one of the long-standing problems 
in dislocation theory. This work represents the first the- 
oretical effort to challenge the traditional point of view 
that regards vacancy as a locking agent for dislocation 
motion. We predict that vacancies can increase partial 
separation distance in Al, and finally we find there ex- 
ists a weak repulsion between dislocations and vacancy, 
which is independent of the dislocation character. 
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FIG. 1: The 7-surface (J/m 2 ) for displacements along a (111) 
plane for (a) pure Al and (b) Al+V systems. The corners of 
the plane and its center correspond to identical equilibrium 
configuration, i.e., the ideal lattice. The two surfaces are dis- 
played in exactly the same perspective and on the same en- 
ergy scale to facilitate comparison. The 7-surface of Al+V is 
truncated to emphasize the more interesting region. 



TABLE I: Fault vectors and energies (J/m 2 ) for some impor- 
tant stacking faults of the pure Al and the Al+V systems. 





Vector 


Al 


Al+H 


Intrinsic stacking 


1/6[121] 


0.164 


0.105 


Unstable stacking 


1/10[121] 


0.224 


0.143 


Unstable stacking 


1/4[101] 


0.250 


0.427 


Run-on stacking 


1/3[121] 


0.400 


0.831 
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FIG. 2: Dislocation Burgers vector density for four disloca- 
tions (clockwise) : screw (0°), 30°, 60° and edge (90°) for the 
pure Al (solid lines) and the Al+V (dashed lines) systems. 
The peaks in the density plot represent partial dislocations. 



TABLE II: Peierls stress (a p , lCT 3 /i), half core width «, 
A), core energies (U CO re, eV/A) for the four dislocations in 
the pure Al and the Al+V systems and binding energy (Ub, 
eV/vacancy) for the four dislocations. 







screw 


30° 


60° 


edge 


dp 


Al 


8.82 


1.77 


3.40 


0.10 




Al+V 


0.69 


0.10 


0.26 


0.05 


c 


Al 


2.1 


2.5 


3.0 


3.5 




Al+V 


3.2 


4.0 


5.6 


6.3 


U C ore 


Al 


-0.084 


-0.110 


-0.168 


-0.198 




Al+V 


-0.046 


-0.073 


-0.133 


-0.164 


Ub 




0.038 


0.037 


0.035 


0.034 



